Crossing the Phantom Divide 
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We consider fluid perturbations close to the "phantom divide" characterised by p = —p and discuss 
the conditions under which divergencies in the perturbations can be avoided. We find that the 
behaviour of the perturbations depends crucially on the prescription for the pressure perturbation 
Sp. The pressure perturbation is usually defined using the dark energy rest-frame, but we show that 
this frame becomes unphysical at the divide. If the pressure perturbation is kept finite in any other 
frame, then the phantom divide can be crossed. Our findings are important for generalised fiuid 
dark energy used in data analysis (since current cosmological data sets indicate that the dark energy 
is characterised by p « — p so that p < —p cannot be excluded) as well as for any models crossing 
the phantom divide, like some modified gravity, coupled dark energy and braneworld models. We 
also illustrate the results by an explicit calculation for the "Quintom" case with two scalar fields. 

PACS numbers: 98.80.-k; 95.36. +x 



I. INTRODUCTION 

The discovery by the supernova surveys P, Q that the 
expansion of the universe is currently accelerating came 
as a great surprise to cosmologists. Within the standard 
cosmological framework of a nearly isotropic and homo- 
geneous universe and an evolution described by General 
Relativity, this behaviour requires a component with a 
negative pressure p < — p/3, commonly dubbed dark en- 
ergy. However, it is very difficult to understand why the 
dark energy should appear at such a low energy scale, 
and why it should start to dominate the overall energy 
density just now. Explaining its nature is correspond- 
ingly regarded as one of the most important problems in 
observational cosmology. 

Current limits on the equation of state parameter 
w = p/ p of the dark energy seem to_indicate that p w —p 
[1,0], sometimes even that p < —p [^], often called phan- 
tom energy [^ . Although there is no problem to consider 
w < —1 for the background evolution, there are appar- 
ent divergencies appearing in the perturbations when a 
model tries to cross the "phantom divide" w = — 1[7]. 
Even though this region may be unphysical at the quan- 
tum level [1, [^ , it is still important to be able to probe 
it, not least to test for alternative theories of gravity or 
higher dimensional models which can give rise to an ef- 
fective phantom energy [l^, [ill, [12, [3 [3- It would 
certainly be unwise to build a strong bias like w > —I 
into our analysis tools as long as experiments do not rule 
it out. In this paper we consider the evolution of the 
perturbations for models where w crosses —1. We find 
that in many realistic cases the divergencies are only ap- 
parent and can be avoided. At the level of cosmological 
first order perturbation theory, there is no fundamental 



limitation that prevents an effective fluid from crossing 
the phantom divide. 

The paper is organised as follows: We start with a 
short recapitulation of the first order perturbations in 
fluids, which also serves to define our notation. In section 
III we study the behaviour of barotropic fluids close to 
p = —p, concluding that this class of fluids cannot cross 
the phantom divide self-consistently. We allow for non- 
adiabatic perturbations in section IV and show that the 
phantom divide can now be crossed as long as we define 
the pressure perturbation in a frame that stays physical. 
We then illustrate the results with the Quintom model of 
two scalar fields before presenting our conclusions. The 
appendices finally discuss the calculation of the effective 
Quintom perturbations in more detail. 



II. FIRST ORDER PERTURBATIONS 

In this paper we use overdots to refer to derivatives 
with respect to conformal time 77 which is related to the 
physical time t hy dt = adrj. We will denote the phys- 
ical Hubble parameter with H and with Ti. the confor- 
mal Hubble parameter. For simplicity, we consider a flat 
universe containing only (cold dark) matter and a dark 
energy fluid, so that the Hubble parameter is given by 



a dt 



-3 



+ (1 - n„,) /(a)] (1) 



implying that the 



where /(a) = cxp —3 J° i^t^^dw 

scale factor today is ao = 1. We will assume that the 
universe is filled with perfect fluids only, so that the en- 
ergy momentum tensor of each component is given by 



{p + p) u'^u^ + p g 



(2) 
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where p and p are the density and the pressure of the 
fluid respectively and u'^ is the four-velocity. This is po- 
tentially a strong assumption for the dark energy, but a 
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more general model should lead to more freedom in the 
dark energy evolution, not less. 

We will consider linear perturbations about a spatially- 
flat background model, defined by the line of element: 



[- (1 + 2 A) dTf + 2B,dr]dx' + 
+ ((1 + 2Hl) % + 2Hnj) dx.dx^] (3) 

where A is the scalar potential; Bi a vector shift; Hl 
is the scalar perturbation to the spatial curvature; 
is the trace-free distortion to the spatial metric, see e.g. 
[m for more details. 

The components of the perturbed energy moment ten- 
sor can be written as: 



^0 — 

j 

rpi 



- (p + dp) 

iP + P) (vj-Bj) 

ip + P) v' 



(4) 
(5) 
(6) 
(7) 



Here p and p are the energy density and pressure of the 
homogeneous and isotropic background universe, 5p is the 
density perturbation, 5p is the pressure perturbation, 
is the vector velocity and 11* is the anisotropic stress 
perturbation tensor [l5|. 

We want to investigate only the scalar modes of the 
perturbation equations. We choose the Newtonian gauge 
(also known as the longitudinal gauge) which is very 
simple for scalar perturbations because they are char- 
acterised by two scalar potentials "0 f^nd 0; the metric 
Eq. ^ becomes: 



ds' 



c? \- (1 + 2V') drf- + (1 - 2(/)) dx.dx"- 



(8) 



and Ul} 



where we have set the shift vector Bi 
0. The advantage of using the Newtonian gauge is that 
the metric tensor g^^i, is diagonal and this simplifies the 
calculations 

The energy-momentum tensor components in the New- 
tonian gauge become: 



^0 



-{p + 5p) 
-ik,T^ ^{p + p)i 
{p + 5p)5] 



(9) 
(10) 

(11) 



where we have defined the variable 6 = ikjV^ which rep- 
resents the divergence of the velocity field and we have 
also assumed that the anisotropic stress vanishes, II j = 
(implying (j) — -tp). 

The perturbation equations are (l6j : 



S = -(l + w) 



— (i-3w)e — 

a 1 



w 



+k 



2 ^P/P 



(12) 



(13) 



As w — > — 1 the terms containing l/{l + w) will generally 
diverge. This can be avoided by replacing 9 with a new 
variable V defined via V — {1 + w) 9. This corresponds 
to rewriting the Oi component of the energy momentum 
tensor as ikjT^ — pV which avoids problems if Tq ^ 
when p — —p. Replacing the time derivatives by a 
derivative with respect to the scale factor a (denoted by 
a prime), we obtain: 



5' 
V 



3(1 



(1 - 3w) 



V 

V P 
a 



1 



5p 



Sp 
P 



wS 



(14) 



Ha^ p 



.)^0.(15) 



In this form everything looks perfectly finite even at w — 
— 1. But in order to close the system, we need to give 
an expression for the pressure perturbations. A priori 
this is a free choice which will describe some physical 
properties of the fluid. In the following we will consider 
several possible choices and discuss how they influence 
the behaviour of the fluid at the phantom divide. We 
will see how the specification of Sp determines if a fluid 
can cross the divide or not. 

In addition to the fluid perturbation equations we need 
to add the equation for the gravitational potential ?/;. If 
there are several fluids present, then the evolution of each 
of them will be governed by their own set of equations 
for their matter variables {Si,Vi}, linked by a common tp 
(which receives contributions from all the fluids) [1^ . We 
use UM. 



^ 2a 



1 fHi^iPi\ ( k 



11*. 



(16) 



The evolution of a fluid is therefore governed by 
Eqs. (fTI)) and (fT5|l. supplemented by a prescription for 
the internal physics (given by dp) and the external 
physics (through -0 and 0'). There are two points worth 
emphasising. Firstly, the dark energy fluid perturbations 
cannot be self-consistently set to zero if w ^ —1. Even 
\i 5 — V = Q on some initial hypersurface, it is unavoid- 
able that perturbations will be generated by the presence 
of ?/;. As this function describes physics external to the 
fluid, it cannot be controlled directly. Conceivably 5p 
could be chosen to cancel the external source of either 5 
or V ^ but not of both. Even then such a 5p would have 
to be incredibly fine-tuned as ip and 0' depend on the 
evolution of the other fluids as well. 

Secondly, we see that as — » — 1 the external sources 
are turned off. A fluid can therefore mimic a cosmological 
constant as then (and only then) 5 = V = 5p = is 
a solution. This corresponds to an energy momentum 
tensor ([2]) with T^^ = pg^v, so that K ^ p. However, in 
general the perturbations do not vanish even if w = — 1. 
A perfect fluid is therefore in general not a cosmological 
constant even if p = —p. 

To calculate perturbations in different gauges we need 
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to introduce the coordinate transformation: 

rj = fj + T (17) 
x' = x' + L\ (18) 

the gauge transformation of the matter variables is then: 



Sp = dp + 3-{l + w)pT 



a 

6p Sp , 



+ 3-n+w)—T 
pa w 



(19) 

(20) 
(21) 



where we have introduced a new quantity — p/ p, cahed 
adiabatic sound speed. 
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III. BAROTROPIC FLUIDS 

We define a fluid to be barotropic if the pressure p de- 
pends strictly only on the energy density p: p = p{p). 
These fluids have only adiabatic perturbations, so that 
they are often called adiabatic. We can write their pres- 
sure as 



P{p) = P{P + 5p) ^ p{p) + 



dp 
dp 



5p + 0{{5pf 



(22) 



Here p{p) = p \s the pressure of the isotropic and ho- 
mogeneous part of the fluid. Introducing A'^ = In a as a 
new time variable, we can rewrite the background energy 
conservation equation, p = —3Ti{p + p) in terms of w, 

dw dp dw dw 
dN^dN^^-^^' + ''^PW 

We see that the rate of change of w slows down as w ^ 
— 1, and w = — 1 is not reached in finite time [13], except 
if dw / dp diverges (or p, but that would lead to a singular 
cosmology). The physical reason is that we demand p to 
be a unique function of p, but at w = — 1 we find that 
p = Q. If the fluid crosses w = —1 the energy density 
p will first decrease and then increase again, while the 
pressure p will monotonically decrease (at least near the 
crossing), Fig. [1] It is therefore impossible to maintain 
a one-to-one relationship between p and p, see Fig. [2] 
(notice that the maximum of p and the minimum of p do 
not coincide). 

An example that can potentially cross w — — 1 is given 

by 



[l + w[p)f = C\-p~-p^). 



(24) 



Starting with w > ~\ and p > px both will decrease 
until w = —1 and p — p^. If it is possible to switch 
from the branch with w > — 1 to the other one at this 
point, then p can start to grow again while w continues 
to decrease. Using the evolution equation ([^5]) for w, we 
find 



dw 
dN 



[C^-px + (1 ■ 



FIG. 1: Energy density (black solid line) and pressure (red 
dashed line) as functions of the scale factor, w crosses —1 at 
ax = 0.5, so that p is minimal at this point while p decreases 
monotonically there (but has a maximum earlier). 



SO that dw/dN ~ — 3/2C^px near the crossing. The full 
solution (with N — at w — —1) is 

wiN) = -l-C^tan(^^^{N-Nx)), (26) 



clearly not a realistic solution for our universe as there 
are divergencies in the finite past and future. Of course 
we could modify the example (j24p . but we are mostly 
interested at the behaviour near the crossing. In our case 
1 -f- w{a) oc —{a — ay ), ie. we observe a linear behaviour. 
This is actually the limiting case: if we choose {1 + w) (x 
{p — pxY close to the crossing, we find 



d{l + w) 
dN 



-C(l + w) 



{2u-l)/v 



(27) 



for small |1 + w\. Correspondingly we cross with zero 
slope for V > 1/2 and with infinite slope for v < 1/2. 
However, if dw/dN = at w = — 1 then the system will 
turn into a cosmological constant and stay there forever. 
On the other hand, this is unstable against small per- 
turbations towards p < —p, so that the crossing might 
eventually be completed anyway. 

Let us have a closer look at the perturbations. The 
second term in the expansion (|22p can be re-written as 



dp 
dp 



P 
P 



3n{l + w) 



(28) 



w 



where we used the equation of state and the conservation 
equation for the dark energy density in the background. 
We notice that the adiabatic sound speed will neces- 
sarily diverge for any fluid where w crosses —1. 

For a perfect barotropic fluid the adiabatic sound speed 
turns out to be the physical propagation speed of 

(25) 

perturbations. It should therefore never be larger than 
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FIG. 2: The pressure as function of the energy density. The 
graph of p{p) shows that p is not a single valued function of 
the energy density, and that there is a point with an infinite 
slope (corresponding to the minimum of p and the divergence 
of c^). The point where p{p) has a zero slope corresponds to 
the maximum of p and a vanishing ci ■ 



The pressure perturbation behaves like 

Sp^p (^(a - ax)"-i + Oiia - ax)")) . (33) 

We need a > 1, otherwise the pressure perturbation will 
diverge as a power-law at the crossing. Unfortunately 
this condition is not sufficient. Looking at the second 
perturbation equation (jlSp we see that although noth- 
ing diverges, there is no reason for V to go to zero at 
the crossing, and in general it will not, except if we fine- 
tune it with infinite precision. As S vanishes at crossing 
and can so potentially cancel the divergence in the sound 
speed, this term is no longer necessarily dominant in the 
differential equation for 6, Eq. (PTjl . We also need to 
take into account the other contributions. The term con- 
taining is sufficiently suppressed by the 1 + w factor 
to neglect it. Taking V to be constant at Cx (to lowest 
order) V{ax) — Vx , we find: 



S' 



-Cx 



(34) 



the sgeed of hght, otherwise our theory becomes acausal 
The condition < 1 for (1 + w) > (our point 



of departure) is equivalent to 



dw 
dN 



<3(l + w)(l-w;). 



(29) 



No barotropic fluid can therefore pass through w = ±1 
from above without violating causality. Even worse, the 
pressure perturbation 



cjp 



w 



(30) 



will necessarily diverge if w crosses —1 and 5p 0. Us- 
ing the gauge transformation Eqs. (fT9| and ([20|l we see 
that the relation between pressure and energy density 
perturbations is gauge invariant. If the pressure pertur- 
bation diverges in one frame, it will diverge in all frames. 
The only possible way out would be to force Sp ^ fast 
enough at the crossing. Let us study the behaviour close 
to i(j = — 1 in some detail now: 

First we look only at the dominant contribution to the 
right hand side of Eq. (fT4|) . Near w = —1 this is clearly: 



S' 



3 2. 

--cj = 
a 



1 



(31) 



Assuming that near the crossing the equation of state 
behaves like 1 4- w « A(a — )" with a > we find the 
solution: 



S (X (a — Ox)" 



(32) 



(independent of A which drops out of the equation). This 
solution goes to zero at the crossing. If the sign had 
been different, we would have found instead the solution 
6 (X {a — Ox)~" which diverges. 



with Cx 
is: 



Ux / (^ (a X ) a x ) ■ The solution to this equation 



c _ / (a - flx) ((5x - Cx fog |a - flx I) a = l . 
1 <5x(a-ax)" + (a-ax)Cx/(l-a) a^l ^^^^ 

We notice that for a > 1 the density perturbations vanish 
at crossing. However, for a = 1 they do not vanish fast 
enough. Instead of the behaviour given in Eq. (|33p. the 
pressure perturbation exhibits now a logarithmic diver- 
gence since the sound speed cancels the factor ( 
Even though our derivation here is not rigorous, numer- 
ical calculations confirm this behaviour, also if the back- 
ground is not completely matter dominated, see Fig. [3l 
Although a full solution may be possible in some cases, 
it would turn out to be a special function, obscuring the 
structure of the result while still showing essentially the 
same simple behaviour. 

We find therefore three possible behaviours for the per- 
turbations, depending on how w crosses the phantom di- 
vide: 

1. a < 1: In this case w crosses —1 with an infi- 
nite slope, leading to a power-law divergence of the 
pressure perturbation. 

2. a — 1: w crosses —1 with a finite, but non-zero 
slope, dp vanishes, V stays finite and generally non- 
zero and dp diverges logarithmically. 

3. a > 1: w crosses —1 with a zero slope. 6p vanishes, 
V and Sp stay finite and generally non-zero. 

The only acceptable case is the last one, but as dis- 
cussed at the beginning of this section, we expect the 
system to get stuck at w — —1 and to never cross. 
Barotropic perfect fluids therefore fail to cross either at 
the background or perturbation level. Even if the fluctu- 
ations can lead to a crossing with zero slope due to the 
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FIG. 3: This figure shows the logarithmic divergence of 
(5/(a — flx) near ax = 0.5. The black curve is the numeri- 
cal solution while the red dashed line shows the approxima- 
tion given by Eq. p5|l for a = 1. The approximate formula 
describes the structure of the divergence quite well. The pres- 
sure perturbation (5p exhibits a very similar divergence in this 
case. 



instability of the cosmological constant solution in some 
cases, the model is still not realistic. The perturbations 
seem to propagate acausally, and generically after the 
transition < for a period, leading to classical insta- 
bilities with exponential growth of the perturbations. In 
the next section we will relax the barotropic assumption, 
which allows for entropy perturbations that can stabilise 
the system and keep the sound speed below the speed of 
Hght. 



the pressure perturbation. However, there are two rel- 
evant variables that we could compare (5p to, the fluid 
velocity V and the perturbation in the energy density bp. 
Clearly it would be counterproductive to replace a single 
free function by two free functions, and it would lead to 
degeneracies between the two. Another problem is that 
the perturbation variables depend on the gauge choice. 
But in this case the two problems cancel each other, lead- 
ing to a simpler solution: Going to the rest-frame of the 
fluid both fixes the gauge and renders the fluid velocity 
physically irrelevant, so that we can now write pol |: 



(36) 



where a hat denotes quantities in the rest-frame. The 
physical interpretation is that (?^(k,t) is the speed with 
which fluctuations in the fluid propagate, ie. the sound 
speed. Again, some physical models lead to simple pre- 
scriptions for the sound speed. The barotropic models 
discussed in the last section have c^{k,t) — c^{t), and 
the perturbations in a scalar fleld correspond (to linear 
order) exactly to those of a fluid with Cg(fc, t) = 1. 

The rest-frame is chosen so that the energy-momentum 
tensor looks diagonal to an observer in this frame. In 
terms of the gauge transformations this amounts to 
choosing B — 0de- However, we notice immediately a 
potential flaw in this prescription close to w — —1: The 
off-diagonal entries of the energy-momentum tensor are 
actually {p+p)9 so that demanding 6—B = is a stronger 
condition than required. In other words, the condition to 
be at rest with respect to the flow of p+p cannot be main- 
tained aX p + p — 0. However, as p > 0, we could define 
instead the sound speed in a frame where there is no fiow 
of energy density. 

To see this, let us calculate the pressure perturbation 
defined by Eq. ([55]) in the conformal Newtonian frame, 
following ^2^: Breaking the single link between p and 
p amounts to the introduction of entropy perturbations. 
A gauge invariant entropy perturbation variable is T = 
f -$f MM- By using 



IV. NON-ADIABATIC FLUIDS 

The discussion of the barotropic fluid shows that we 
have to violate the constraint that p be a function of p 
alone. At the level of first order perturbation theory, 
this amounts to changing the prescription for Sp which 
now becomes an arbitrary function of k and t. This 
problem is conceptually similar to choosing the back- 
ground pressure p(t) , where the conventional solution is 
to compare the pressure with the energy density by set- 
ting p(t) = w{t)p(t). In this way we avoid having to 
deal with a dimensionfuU quantity and can instead set 
w, which has no units (up to a factor of the speed of 
light) and so is generically of order unity and often has a 
simple form, like w = 1/3 for a radiation fluid or w = —1 
for a cosmological constant. 

It certainly makes sense to try a similar approach for 



wT 



p f Sp Sp 



(37) 



and the expression for the gauge transformation of Sp [2]| , 



_V 

Sp = Sp + iHp-pr 



we find that the pressure perturbation is given by 



Sp = cl 



V 



(38) 



(39) 



As i 



oo at the crossing, it is impossible that all other 



variables stay finite except if V ^ fast enough. Again 
we will show that this is not in general the case, except if 
— > or ^ at crossing. However, we will then see 
that this is not required, indeed we will argue in the next 
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section that the more generic sohition is to let diverge 
at the crossing in order to cancel the divergence of c^l 

To this end, we consider again the structure of the 
perturbation equations near w = —1. Inserting the ex- 
pression for dp into our system of perturbation equations 
we find 



S' = 3(l + u;)V/ 
3 



V_ 

fc2 



fc2 



+■ 



w 



(40) 



(41) 



The spoiler here is the continued presence of c„ which we 
know to diverge at crossing. Let us start by considering a 
finite and constant . For all reasonable choices of k the 
dominant term in the V equation will be the one contain- 
ing c^. We find that this time the velocity perturbation 
is driven to zero at crossing: 



v 



1 + w 



-V. 



(42) 



Proceeding as in the last section, we find that now V cx 
(a — flx)", cancelling the divergence in for a> 1. Now 
Sp will in general not vanish at crossing, and we have to 
include that term in the differential equation for V . The 
term with ijj is suppressed by 1 + w and is of higher order. 
As in the last section for S, the solutions for V are now 
to lowest order: 



V 



(a - ax) (Vx - -Dx log |a - ax I) a = l 
l^x(a-flx)" + (a-ax)i5x/(l-a) a^l 



(43) 

with Dx being clk^d/{Ha^), evaluated at crossing. 

Again there are no divergences appearing in the energy 
momentum tensor only if a > 1, ie. if u; crosses —1 with a 
zero slope. A possible way to get around this fine-tuning 
is to demand that = at crossing, as in this case 
the logarithmically divergent term disappears. We also 
notice that the usual velocity perturbation 9 = V/(l + w) 
does diverge in all cases, either logarithmically if a = 1 or 
as (a — ax)^~" if a > 1. For an observer in the rest-frame 
where B = 9 this means that the metric perturbations 
become large - at crossing the metric is even singular. 
At the very least, perturbation theory is no longer valid 
for such an observer. 

Another way to see that the rest-frame is ill-defined is 
to look at the energy momentum tensor For p — —p 
the first term disappears, leaving us with T'^'^ — p g^''. 
Normally the four- velocity is the time-like eigenvector 
of the energy momentum tensor, but now suddenly all 
vectors are eigenvectors. The problem of fixing a unique 
rest-frame is therefore no longer well-posed. 

However, hy construction the pressure perturbation 
looks perfectly fine for precisely the observer in the rest- 
frame, as 5p = c^(5/5 does not diverge. Our prescription 



for the pressure perturbations has singled out the one 
frame which we cannot use for fluids crossing the phan- 
tom divide. The reason is that the gauge transformation 
relating the pressure perturbations in the different gauges 
is: 



5p = 5p + SHpc^^ 



V_ 

fc2' 



(44) 



If V does not vanish fast enough at the crossing then 
the pressure perturbation has to diverge in at least one 
frame. As we have just discussed, the dark-energy rest- 
frame becomes unphysical at crossing. It is clearly better 
to specify a finite pressure perturbation in a different 
frame. 

One problem is to find a way of characterising the pres- 
sure perturbations in a physical way - in the Quintom 
example of the next section, we find that the additional 
contributions diverge but in such a way that we end up 
with a finite result for dp. As another example, we just 
choose dp proportional to 5p, 



5p{k, i) = 7 5p{k, t) 



(45) 



in the conformal Newtonian gauge. This will work as 
long as i5p ^ 0, otherwise it forces 5p to vanish in the 
same place as 5p, which is not general enough. If we 
insert this expression into eqs. (fTil) and (fT5|) . we obtain: 



5' = 3(1 + 1^)7/.' --^-3^(7-^)5 (46) 

a Ha^ Ha"' 

As u; ^ — 1 none of the terms diverge, so that S and V 
stay in general finite and non-zero at crossing. We show 
in Figs. |l]and[S]a numerical example for two choices of 
7 where it is impossible to see that the phantom divide 
has been crossed at ax = 0.5. 

It is of course possible to express 7 in terms of cj., and 
in terms of 7. Using the expression for the pressure 
perturbation in the rest-frame given by Eq. (j36p and the 
gauge transformation given by Eqs. ([TO]) and ([^0]) . we 
obtain 



7(5p + SHcIp^ 



Sp 



V 



(48) 



In general 7 will therefore be scale dependent even if 
cl is not (even though of course will also in general 
depend on k and t), or vice versa. Also, to reproduce the 
evolution with finite 7 shown in Figs. [3] and we would 
have to substitute a divergent to cancel the divergence 
in (cf. figure [7| which shows how the apparent sound 
speed diverges in the Quintom example). We also notice 
that on very small scales where fc ^ 1 we find 7 ~ for 
finite c^, which is the usual result that on small scales 
gauge differences become irrelevant (for physical gauges). 

Finally we would like to emphasise again that for a 
general, finite pressure perturbation in any gauge except 
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FIG. 4: This figure shows the density contrast Sde (red and 
blue lines, negative for small scale factor a) and the velocity 
perturbation Vde (cyan and magenta lines, positive for small 
a) for a dark energy component with a pressure perturbation 
given in the conformal Newtonian gauge through 7 (as defined 
in Eq. psp ). The solid lines show the results for 7 = 0.2 and 
the dashed lines for 7 = l.The crossing of the phantom divide 
at flx = 0.5 (dotted green vertical line) is not apparent in 
these figures, everything is finite. 



FIG. 5: The gravitational potential 7/) in the same scenario as 
Fig. |4] (with 7 = 0.2 for the black solid line and 7 = 1 for the 
red dashed line). Again, the crossing of the phantom divide 
at flx = 0.5 (dotted green vertical line) is not apparent. The 
gravitational potential is constant during matter domination 
(small a) and starts to decay as the dark energy begins to 
dominate and the expansion rate of the universe accelerates. 
The contribution of the dark energy perturbations to tp is very 
small {<,0.1%) so that there is no visible difference for the 
two choices of 7. 



the fluid rest-frame, there is no problem with the per- 
tm'bation evolution across the phantom divide. Also, in 
general the perturbations, including V, will not vanish 
at crossing. 



V. THE QUINTOM MODEL AS EXPLICIT 
EXAMPLE 

To clarify some of the above points we consider an 
explicit example of a model crossing the phantom divide, 
the Quintom model p^ . l23j , and compute the pressure 
perturbation ab initio at crossing (see also [24|). 

The original Quintom model considered two scalar 
fields. For us it is advantageous to use instead two flu- 
ids with a constant and equal rest-frame sound speed 
Cg = 1. At the level of first order perturbation theory 
the two models are exactly equivalent, see appendix [Xl 
As in the original model, we use two fluids with constant 
equations of state parameters wi and 'W2- We will have 
to map the two fluids onto a single effective fluid. To this 
end we define the effective parameters in such a way that 
the effective energy momentum tensor is the sum of the 



two fluid energy momentum tensors. This leads to: 



Pcff 


= Pl+ P2 






(49) 


Pcff 


= PI+P2 






(50) 


WcS 


Wipi + W2P2 
Pi + P2 






(51) 


5cfi 


PlSl + P2S2 
Pi + P2 






(52) 




(1 + wi) pi9i - 


Ki- 


1- W2) P2^2 


(53) 


(1 -1- Wi) Pi - 




VW2)P2 



Re-expressing the perturbation equations in these vari- 
ables we find precisely the normal perturbation equations 
for a single fluid and p5|) with pefi and Jpeff re- 
placing p and 5p. PcB is simply given by WeSPeS, but 
the density perturbation is more involved. Starting from 
SpcS = Spi + dp2 we can write it as 

6pcS = Cs off ^Pcff + Sprc\ + Spnad 

+3n{cl,g-cl)p,s^ (54) 

where the (effective) adiabatic sound speed is as always 
= pcs/pcS and where ^ff — 1 for two scalar fields. 
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The other terms are the relative pressure perturbation 

(55) 



{W2 - Wi) piP2 

oPici = 7::rr^ "^12 



SHpcff 

given by the relative density perturbation of the two 
scalar fields, 



5*12 — 



Si 



1 + Wi 1 +W2 



(56) 



(which corresponds to a gauge invariant relative entropy 
perturbation [25j) as well as the non-adiabatic term 



Spnad = - {W2 - Wi) 



P1P2 



S 



12 



(57) 



given by the relative motion of the two scalar fields with 
A9i2 = 9i — 62 (see appendix ICj). 

We notice that the relative perturbations act as inter- 
nal perturbations which couple to the effective variables 
purely through the behaviour of the pressure perturba- 
tion. Understanding the physics behind the dark energy 
in such a case will therefore require a precise measure- 
ment of the pressure perturbations as well as their careful 
analysis, to uncover the different internal contributions. 

Even though the effective sound speed = 1 in 

Eq. ([54|) is finite in the rest-frame, the transformation to 
any other frame will lead to a divergence as discussed in 
the previous section. This divergence then needs to be 
cancelled by the two additional terms, (5proi and (5pnad- 
In our example, both diverge, cancelling together the di- 
vergent contribution from the singular gauge transforma- 
tion as well as their own divergencies, see Fig. [51 This 
behaviour, which looks extremely fine-tuned, is automat- 
ically enforced in this model. Such a cancellation mecha- 
nism is required for any model in order to cross w = —I. 

We also see that although the effective sound speed 
remains simply = 1, this is only true if we know that 
there are internal relative and non-adiabatic pressure per- 
turbations (as well as their form). But in general we 
would try to parametrise the pressure perturbation as: 



6pcs = clSpcs + 3H (cl 



Vr, 



PeS^. (58) 



where the apparent sound speed is now a mixture of 
the real effective sound speed together with the relative 
and the non-adiabatic pressure perturbations. In this 
case we no longer find a simple form for the sound speed. 
In Fig. [7] we plot as a function of the equation of state 
parameter w for several wave vectors k. As predicted, the 
apparent rest-frame sound speed diverges at the crossing. 
In the Quintom case the effective perturbations do not 
vanish at w = — 1. 

We think that the lessons learned from the Quin- 
tom model are applicable also to more general models 
with multiple fields, non-minimally coupled scalar fields, 
brane-world models and other modified gravity models 
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FIG. 6: This figure shows the different divergent contributions 
to the pressure perturbation, Eq. (|54|) . multiplied by lO''. The 
relative pressure perturbation is shown as red dotted line, the 
non-adiabatic pressure perturbation as green dash-dotted line 
and the contribution from the gauge transformation to the 
conformal Newtonian frame as blue dashed line. Each of the 
contributions diverges at the phantom crossing, but their sum 
(shown as black solid line), and so 5p, stays finite. 




FIG. 7: We plot the apparent sound speed defined by 
Eq. (|58p for three different wave vectors, k = 1/Ho (black 
solid line), k = 10/ Hq (red dashed line) and lOO/ffo (blue 
dotted fine). Although the real sound speed is just = 1, 



the apparent sound speed diverges at w = 
become negative. 



-1 and can even 
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that can be represented by an effective dark energy fluid. 
In all these models, as in the Quintom case, WcS = ~1 
is not a special point in their evolution. There is no rea- 
son to expect that the model will adjust its behaviour at 
this point, as the crossing of the phantom divide is in- 
cidential. In representing these kinds of models with an 
effective fluid we therefore expect that the perturbations 
will not vanish, and that the apparent sound speed would 
have to diverge. 

One important difference between the Quintom model 
and more general modified gravity models is that the 
latter seem to require generically a non-zero anisotropic 
stress as if) ^ ip [2^. In this case — does not allow 
crossing the phantom divide, since in Eq. ([43]) con- 
tains an additional non-zero term due to the anisotropic 
stress. This reinforces our view that the Quintom-like 
crossing, where all perturbations stay non-zero and SpoS 
remains finite in spite of divergences in most convention- 
ally expected terms, is more generic. This means that 
there is no obvious way to predict the form of SpeS in gen- 
eral for modified gravity models. On the positive side, if 
the dark energy is not just a cosmological constant, and if 
we are able to measure SpcS then we may hope that this 
will provide us with clues about the physical mechanism 
that is causing the accelerated expansion of the universe. 



VI. CONCLUSIONS 

As long as the cosmological data indicates the pres- 
ence of a dark energy with an effective equation of state 
p « —p it will be necessary to consider models with the 
same equation of state. In general, we will have to al- 
low the equation of state to cross the phantom divide, 
w = —I. Even though such a fluid model may not be 
viable at the quantum level, it is possible that this be- 
haviour is only apparent, or due to a modification of 
General Relativity or the existence of more spatial di- 
mensions. In analysing the data we therefore have to be 
able to use general self-consistent models at the level of 
linear perturbation theory. 

In this paper we have studied the behaviour of the per- 
turbations in general perfect fluid models close to w ~ 
— 1. We have shown that although models with purely 
adiabatic perturbations cannot cross w = —1 without 
violating important physical constraints (like causality 
or smallness of the perturbations), it is possible to rec- 
tify the situation by allowing for non-adiabatic sources of 
pressure perturbations. However, the parametrisation of 
dp in terms of the rest-frame perturbations of the energy 
density cannot be used as this frame becomes unphysical 
at w = —1. By parameterising Sp instead in any other 
frame the divergencies are avoided. 

We also computed all quantities in the Quintom model 
which provides an explicit example of the above mecha- 
nism. In this model, even though the propagation speed 
of sound waves remains finite and constant, the addi- 
tional internal and relative pressure perturbations lead 



to an apparent sound speed which diverges. It is only 
the sum of all contributions to dp which remains finite. 

A more speculative conclusion is that it seems difficult 
for "fundamental" fields to cross w = — 1 as their ap- 
parent rest- frame sound speed, defined through Eq. (j36p . 
would generally be the actual propagation speed of their 
perturbations, which must remain smaller than the speed 
of light. From the discussion in section [TVl we learn that 
fields with a well defined rest-frame and sound speed have 
only two routes to phantom crossing: 

• dw / da = at crossing: This looks rather fine-tuned 
as the field needs to be aware of the presence of the 
barrier. On the other hand, a normal minimally 
coupled scalar field reaches w — ~1 when </> = 0, 
and it does so with zero slope. A scenario where 
phantom crossing of this kind is realised could be 
built by changing the sign in front of the kinetic 
term whenever (j) = 0. As an example, a cosine 
potential then leads to a "phaxion" scenario, see 
Fig. [8l On average, the phaxion behaves like a 
cosmological constant, but oscillates around w = 
— 1. However, it is unclear how to construct such a 
scenario in a covariant way. 

• vanishing sound speed: A field with a vanishing 
sound speed can also avoid the logarithmic diver- 
gence of 6p. This requires either a coupling between 
l + w and or else the sound speed might be zero 
at all times, see Fig. [S) In the latter case the pres- 
sure of the field would need to remain close to —p at 
all times in order to prevent the small-scale pertur- 
bations from growing too quickly. Also, the sound 
speeds needs to be exactly zero, otherwise the log- 
arithmic divergence reappears. This may require a 
symmetry enforcing = as even a very small but 
non-zero sound speed would render this scenario in- 
viable. 

In view of these difficulties, the Quintom family of mod- 
els, although more complicated than a single field, may 
prove to be the conceptually simplest way to cross w = 
— 1 with a model defined through an action. They also il- 
lustrate that the effective pressure perturbation may pro- 
vide a kind of fingerprint of the mechanism behind the 
accelerated expansion of the universe if it can be mea- 
sured. 

Concerning the use of phantom-crossing fluids for the 
purpose of data analysis, it is straightforward to imple- 
ment them by just avoiding the usual parametrisation 
of the pressure perturbations in terms of the rest-frame 
sound speed. However, there does not seem to be a 
canonical way to choose the pressure perturbations if one 
is not allowed to use the dark energy rest-frame. The aim 
for the far future will be to directly measure the pres- 
sure perturbations of the dark energy in order to gain 
insight into the physical origin of the phenomenon. For 
now, we have to ensure that the definition of Sp does not 
lead to unphysical situations, while preserving the usual 
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FIG. 8: The phaxion: The sign in front of the kinetic term of a 
normal scalar field is fiipped every time (j) passes through zero 
(so that > 0). As potential we use V{<j)) = 1 + cos(0/mo) 
which is bounded both from above and below. The field 
then moves to large values (top left graph), while w oscillates 
around —1 (top right). The energy density p also oscillates 
so that the field mimics an effective cosmological constant 
(bottom left). As the phantom divide is crossed with zero 
slope, there is no divergence in the pressure perturbation (bot- 
tom right), but the time-dependent effective mass can lead to 
strong growth of the perturbations in spite of = 1. 



parametrisation in terms of the rest-frame sound speed 
as much as possible far away from the phantom divide. 
Maybe the simplest way out is to regularise the adiabatic 
sound speed which appears in Eq. ([55]) because of the 
gauge transformation into the rest-frame. While any fi- 
nite choice of Sp is a physically acceptable choice, it is 
preferable to modify Eq. (|39p in a minimal way so that 
the usual interpretation of is preserved for w —1. 
We propose to use 



w{l 



3H[(l + w)2 +A] 



(59) 



where A is a tuneable parameter which determines how 
close to ui = —1 the regularisation kicks in. A value 
of A w 1/1000 should work reasonably well, as shown 
in Fig. [ini In this case Sp is well-defined and there are 
no divergencies appearing in the perturbation equations 
(|40| and (|4T|l . Although the differences in Sp for the dif- 
ferent choices of A in Fig. (TU] look important, we have to 
remember that, firstly, the perturbations in the dark en- 
ergy are normally small (especially close to the phantom 
divide) and that secondly they are only communicated 
to the other fluids via the gravitational potential ip- As 
in Fig. [HI we find also here that the dark energy pertur- 
bations are subdominant. Computing the CMB power 



FIG. 9: The pressure perturbation Sp for three fluids with 
different rest- frame sound speed, — 0.1 (top curve, red 
dash-dotted line), = 0.01 (middle curve, blue dashed line) 
and = (lowest, solid line). Only the last case does not 
have a logarithmic divergence in the pressure perturbation at 
the phantom divide. 



spectrum for A = 1/100 and A = 10^^ we find a relative 
difference in the Cg of less than 1% on large scales, which 
rapidly drops to 10~^ by = 50, well below cosmic vari- 
ance. We expect detectable differences only for strongly 
clustering dark energy with a sound speed that is close 
to zero or even negative. Although it is reassuring that 
directly measureable quantities are not sensitive to the 
precise value of A in Eq. ([5^ , this also shows that it will 
be very difficult to distinguish between different models 
of dark energy if w w — 1 . 
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APPENDIX A: EQUIVALENCE BETWEEN 
SCALAR FIELDS AND FLUID MODELS 

The aim of this appendix is to show that at the level 
of first-order perturbation theory a scalar field behaves 
just like a non-adiabatic fluid with = 1. To this end 
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FIG. 10: The pressure perturbation for a fluid example with 
Cs = 0.1. The red dash-dotted curve shows the standard 
case with the logarithmic divergence at w — —1. The blue 
(dashed), black (solid) and green (dotted) curve use a regu- 
larised adiabatic sound speed which does not diverge. The 
black curve with A = 1/1000 provides a reasonable flt. Larger 
values smooth too much, while smaller values start to follow 
the divergence and exhibit a temporary instability in the so- 
lution for V. 



we decompose the scalar field into a homogeneous mode 
4>(t) and a perturbation S(j){k, t). At the background level 
we find 



1 

1 
2a 



4>^ + v [cf) 



(Al) 
(A2) 



and the equation of conservation is just the equation of 
motion, 



2H(l) + a^— ^ 0. 



(A3) 



The adiabatic sound speed is defined as 
3aH 12 _ 2^ ' 

P 



(and we remind the reader that a/a = Ti = aH). The 
perturbed energy momentum tensor is: 



ST! 



-ikST^ = ikST° 







dV 








- ^02^ - 


dV 




" d0 



fc2 
^2 



pV 



(A7) 



where we wrote ^' for the gravitational potential ip in 
order to avoid confusions with the scalar field variables 
(only in this appendix). 

In order to derive the rest frame sound speed of the 
scalar field we use equation ((39|) 



(A8) 



Sp ^ clSp +^{cl- cl) pV 



and express everything in terms of scalar field quantities. 
We find 



02^_ 

dm 



1 



dV . , ZaH 



iaH ( 2a 1 dV\ 



(A9) 



which after some algebraic manipulations turns into 



1 



1 



dV ^ 3aH 



d(f) 



= cl —cPS(f> - —02* + —6cP + -^cPS(f> I 
V d(p / 



AlO) 



Therefore = 1. 

Now let us derive the equation of motion for the scalar 
field perturbations from the perturbation equation for a 
perfect fluid, Eq. 

S = -{l + w)(9-3i'^-3aH(^^~wS^ (All) 

which can be rewritten as 

Sp + 3aH {j5 + p)S = 3{p + p)i' - pV (A12) 

Expressing the time derivative Sp in terms of scalar 
field quantities, 

2ai/ • • 1 •• • 2aH ■■ dV ■ 

Sp = —<t>5^ + —^^4' - —yr<l>^<P - -ttS(I> - (A13) 



a' 



a^ 



a^ 



1 6aH -^^ ^-dV ^ dV ■ d^V 

d(p d(p dcp'^ 



and doing likewise with the other terms, 



iaH {Sp + Sp) 



a 

3 • ■ 

3{p + p)^ = ^0^* 
a-' 

pV = ^ ^ 



6aH ■ 6aH -o / , -, ,n 
^-r^' (A14 



(A15) 
(A16) 



we can insert all these expressions into Eq. (|A12p and 
obtain finally 



2aHS(j) + 



d^V fc2 



dV 



dcj)^ a2 

(A17) 

which is indeed the equation of motion for S(f> in the con- 
formal Newtonian gauge (see e.g. [lB|). 
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APPENDIX B: EFFECTIVE PERTURBATIONS 
IN TWO BAROTROPIC FLUIDS 

We consider first two barotropic perfect fluids with 
constant equation of state parameters, 



pi = wipi, wi > -1; 

P2 = W2P2, W2 < -1. 



(Bl) 
(B2) 



We can define the effective quantities as those appearing 
in the sum of the two energy-momentum tensors, 



PcS 


= Pl+ P2 






(B3) 


PcS 


= PI+P2 






(B4) 




WlPl + W2P2 
Pi +P2 






(B5) 




PlSl + P2S2 
Pi + P2 






(B6) 


ecff 


{l + wi)pi9i^ 




h W2) P2O2 


(B7) 


(1 + wi) pi H 




\- W2) P2 



The system above is characterized by four variables ^i, 
^21 ^ii ^2- In order to have a complete mapping, we 
need to introduce two more variables which express two 
internal degrees of freedom [2^ : 



'5'l2 

M12 



Si 



1 + Wi 1 + W2 
61 — 62 



(B8) 
(B9) 



called respectively, relative entropy perturbation and the 
relative velocity of the two fluids. 

We can now calculate the perturbation equation for 
the effective fluid from 



rpfj, 



(BIO) 



The Tq.^ and T^^^ components give the perturbation 
equations: 

PiSi + P2S2 - 30 [(1 + wi) pi + {1 + W2) P2] + 

+ [(1 + wi) piOi + (1 + W2) P2O2] + 

+2>li [5pi — wipiSi + Sp2 — W2P2S2] = (Bll) 

the derivative of ScS is: 

■ P1S1+P2S2 ^^wipiSi + W2P2S2 , 
Ocff = — : — on- 



Pi + P2 Pl+ P2 

wipi + W2P2 Pi5i + P2S2 



-m 



Pl + P2 



Pi + P2 



(B12) 



inserting the last one into Eq. (jBlip and remembering 
Eqs. (Hil to (|B8|, we have: 



(1 + Wcs) [d^s - 3-0 
5pcs 



Poff 



(B13) 



where 5pcti = Spi + 5p2- This is just the equation (|12|) 
for the effective quantities. 

For the second perturbation equation we have: 



(1 + wi) P161 + (1 + W2) P2Q2 + 
+4W [(1 + wi) piBi + (1 + W2) P2O2] + 

\{i + wif pidi + (1 + W2f P2O2 



-k^lP [(1 + Wi) Pl + (1 + Wi)pi] + 

-e {Spi + 5p2) = 
the derivative of 9cs is: 

(1 + Wi)pi9i + (1 + W2) P2(- 



(B14) 



^cff 



(1 + Wl) Pl + {1+ Wl) Pl 



+ 



-m 



(l + wi) pi9i^{\ + W2) P262 

(1 +Wl)pi + (1 +Wl)pi 
nnj PlOl + (1 + W2) P2d2 

-on — — — — X 

(1 + Wl) Pl + (1 + Wl) Pl 

^ (1 + Wi)^ Pl + (1 + W2f P2 
(1 + Wl) Pl + (1 + Wl) Pl 



(B15) 



again inserting the last one into Eq. (|B14[) and remember- 
ing Eqs. (jB3p to (jB7p and (|B8p . we obtain the expression 
for the second perturbation equation: 



7cff 



(1 - 3wcff) - 
5pcs/pc« , 



+k 



1 



Wcff 



= 



(B16) 



Again, this is the same equation as the one for a single 
perfect fluid, Eq. (fT3)) . The difference to the general case 
is that now the pressure perturbations are fixed by the 
barotropic nature of the two fluids. Starting from the 
generic expression 

6paS = 5pi + 5p2 ^ WipiSi + W2P252 (B17) 

we find that the effective pressure perturbation is 



5pcS = c^<5off - 



1 (W2-Wi)pip2 



512 



(B18) 



37Y pofi 

where is the adiabatic sound speed and is given by 
wipi -I- W2P2 



PcS 

PcS Pl + P2 

Wl (1 + Wl) Pl + W2 (1 + W2) P2 

(1 -t- Wi)pi -I- (1 -I- W2)p2 



(B19) 



The second term appearing in Eq. (jBlSp can be con- 
sidered as the relative pressure perturbation due to the 
relative motion: 



SPtcI = 



(W2 - Wl) P1P2 

STipcS 



(B20) 
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A single barotropic fluid has the pressure perturbation 
Sp — c^Sp. For the two-fluid case we find an additional 
part coming from the relative perturbations of the two 
fluids. The new variable 5*12 given by Eq. (IBSP is a gauge 
invariant relative entropy perturbation [25] . 

The time evolution of S12 is given simply by: 

5i2 = -A012, (B21) 

and contains the fourth variable A9i2, which is a relative 
velocity perturbation and evolves according to 

dtAei2 = -HAOu + {wi + W2 - cl) A012 + 

+ {Wi +W2- cl) Si2 + 

+m K-^;2)geff + fc' r~'"' (B22) 

1 + Wcff 



What we need to evaluate again in the Quintom model 
is the effective pressure perturbation, because all the 
other terms are the same. Taking the sum of the pressure 
perturbation defined in the rest-frame, we have: 

SpeS = C^s,l^Pl + Cs,2'5P2 + 

+mil+W2){ci2^W2)p2^ (C5) 

inserting the Eqs. - (IC4| in Eq. (IC5| we find: 

SpcS = clcffSpcS + Sprd + Spnad 

+3H {cl,g - cl) (1 + w,s) Poff ^ (C6) 



APPENDIX C: EFFECTIVE PERTURBATIONS 
IN THE QUINTOM MODEL 

Perturbations in barotropic fluids with constant w < 
grow very rapidly due to the imaginary sound speed, 
= w < 0. A realistic model crossing the phantom 
divide needs therefore to be composed of fluids with non- 
adiabatic fluctuations and a positive sound speed. In the 
case of the Quintom model, we are dealing with two fluids 
with cl ~ 1. As in appendix [Bl we define the effective 
quantities via eq. (|B3|) - eq. (|B9|) : using these equations 
we find the relations between the variables of the two 
fluids and the variables of the effective fluid: 



61 



62 = 



Soff 



l + Wi 
1 -I- WcS 

(1 -I- Wl) (1 -I- W2) (Wcff - Wl) 
(1 -I- Wcs) iw2 - Wl) 



S 



12 



'JcS 



1 + W2 
1 + WeS 

(1 -I- Wl) (1 + W2) (W2 - Wcff) 
(1 -I- Woff) iw2 - Wl) 

{1 + W2) {w^s - Wl) 

fcff + 7^— TT r^t7i2 



72 = 



(1 + Woff) (W2 - Wl) 
(l+Wi) (Wcff ~ W2) 
(1 -I- Wcff) {W2 - Wl) 



AO 



12 



(Cl) 

(C2) 
(C3) 
(C4) 



where is the effective rest-frame sound speed; Jprci 
is the relative pressure perturbation and (5pnad is the non 
adiabatic contribution to the pressure perturbation; they 
are: 



•-sxfr 

5Prc\ 
<5Pnad 



Cll (1 + Wl) pi + Cj ^ (1 + W2) P2 

{l + Wi)pi + {1 + W2)P2 
{W2 - Wl) P1P2 

SHpcs 



S 



12 



-^s,2j 



{W2 - Wl)] 



P1P2 
SHpeS 



Si2 + -^A9i2 



(C7) 
(C8) 

(C9) 



The total number of degrees of freedom remains the 
same when we change from the "two fluid" to the "single 
effective fluid" picture. In both cases we have four vari- 
ables. What changes is the way these variables interact. 
In the two fluid case, the interaction proceeds through 
the gravitational potential ip- In the single effective fluid 
picture, the additional degrees of freedom become inter- 
nal and appear through additional contributions to the 
pressure perturbation 5p. 
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